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Abstract. A possibly fruitful extension of conventional random matrix ensembles 
is proposed by imposing symmetry constraints on conventional Hermitian matrices or 
parity-time- {VT-) symmetric matrices. To illustrate the main idea, we first study 
2x2 complex Hermitian matrix ensembles with 0(2) invariant constraints, yielding 
novel level-spacing statistics such as singular distributions, half-Gaussian distribution, 
distributions interpolating between GOE (Gaussian Orthogonal Ensemble) distribution 
and half Gaussian distributions, as well as gapped-GOE distribution. Such a 
symmetry-reduction strategy is then used to explore 2x2 PT-symmetric matrix 
ensembles with real eigenvalues. In particular, T'T-symmetric random matrix 
ensembles with U{2) invariance can be constructed, with the conventional complex 
Hermitian random matrix ensemble being a special case. In two examples of VT- 
symmetric random matrix ensembles, the level-spacing distributions are found to be the 
standard GUE (Gaussian Unitary Ensemble) statistics or "truncated-GUE" statistics. 
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1. Introduction 

Random matrix theory (RMT) has found enormous applications in a variety of research 
areas lU |2l |3l H]. Predictions of RMT are often widely applicable and it is of 
fundamental interest to explore its variants and extensions. Two best-known random 
matrix ensembles are Gaussian Unitary Ensemble (GUE) of complex Hermitian matrices 
and Gaussian Orthogonal Ensemble (GOE) of real symmetric matrices. Either GUE 
or GOE can be characterized by P{H)diJ,{H), where yu(i?) is the measure of random 
matrices H and P{H) is the associated probability density. P{H)d^{H) is assumed 
to be invariant under certain symmetry groups, called the symmetry invariant group 
below. The existence of a symmetry invariant group is an important element of RMT, 
indicating the representation-independence of an ensemble within its symmetry class. 
That GUE and GOE represents two different universality classes of statistics is exactly 
because their symmetry invariant groups are different. As seen below, the concept of 
symmetry-invariant group becomes even more important as we extend the conventional 
random matrix ensemble from Hermitian matrices to parity-time- {VT-) symmetric 
matrices. 

One can also define "canonical groups" of transformations for complex Hermitian 
matrices and for real-symmetric matrices f2]. Adopting the terminology in Ref. [2], 
a canonical group for a complex Hermitian matrix is defined as the largest group of 
transformations that can maintain the matrix eigenvalues as well as the Hermitian or 
real-symmetric form. For N x N complex Hermitian matrices, the canonical group is 
U{N)] whereas for N x N real- symmetric matrices, the canonical group is 0{N). In 
general, it is often implicitly assumed that the symmetry invariant group of P{H)dn{H) 
must be identical with the canonical group. This assumption does not always reflect 
a true physical situation and hence can be lifted. For example, if a complex system 
is subject to a known external influence which breaks the maximal symmetry of the 
ensemble, then it is more appropriate to model the system by an ensemble with 
reduced invariant group, e.g., on a subgroup level. Indeed, in earlier studies of the 
GUE-GOE transition of random complex Hermitian matrices [5|, each matrix H is 
decomposed into a real-symmetric matrix plus an anti-symmetric matrix; then based on 
this decomposition a probability density PGUE-GOE(-f^) is proposed, which is invariant 
only under an orthogonal group, not a unitary group. 

The well-known probability density PGUE-GOE(-f^) for GUE-GOE transitions is just 
one possible measure to model an ensemble of complex Hermitian random matrices 
with reduced symmetry. Because 2x2 Hermitian random matrices are just one special 
case of 2 X 2 "PT-symmetric random matrices [6], it should be interesting to explore 
novel random matrix ensembles by proposing non-conventional P{H)dfi{H) with a 
reduced symmetry invariant group and then extend our approach to "PT-symmetric 
random matrix ensembles. Our basic procedure is in fact intuitive: Reducing the 
invariance symmetry of P{H)dn{H) can be realized by imposing a certain constraint 
that is invariant under a reduced symmetry. Though our explicit results are based 
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on 2 X 2 matrices, they can be already relevant to certain physical situations, 
e.g., they may capture short-ranged statistical properties of large random matrices. 
Another motivation to work on 2 x 2 matrices is that the level-spacing distribution 
functions, which reflects the level repulsion properties of random matrices, can be 
obtained analytically for 2x2 matrices with elementary tools only. For random 
Hermitian matrices with reduced invariance symmetry, various interesting results 
of level-spacing statistics are found, including singular distributions, half-Gaussian 
distribution, "gapped-GOE" distribution, and distributions interpolating between GOE 
and half-Gaussian distributions. 

Using our approach to the treatment of Hermitian random matrix ensembles 
with reduced invariance symmetry, we move on to construct ensembles of 2 x 2 VT- 
symmetric random matrices. A general statistical measure is deduced from a statistical- 
independence requirement and then various constraints can be imposed to obtain 
particular cases with reduced invariance symmetry, yielding interesting level statistics 
for non-Hermitian random matrix ensembles (such as a "truncated-GUE" distribution). 
Two important results are as follows. First, the familiar GUE of 2 x 2 random matrices 
can be regarded as a symmetry-reduction result from 2x2 PT-symmetric random 
matrices. Second, there exists a class of 2 x 2 non-Hermitian (but PT-symmetric) 
random matrix ensembles whose level-spacing distribution is still given by GUE. 

Before ending this section, we wish to point out that extending RMT beyond 
Hermitian matrix ensembles has been of interest for decades [71 [2] . Thanks to the newly 
established topic of "PT-symmetric quantum mechanics, this work offers an unforeseen 
route of extension by exploiting a recent study that gives fully-parameterized 2x2 
complex matrices with real eigenvalues [6] . It is hoped that more interests in a broader 
class of random matrix ensembles can be motivated. 



2. General consideration for Hermitian matrices 

A 2 X 2 Hermitian matrix H has four real parameters, 

^ _ ( ^11 "12 + i^l2 

The natural choice of measure dfi{H) = daiidai2da22d6i2 is invariant under the 
canonical group U{2). Under the same symmetry, there are two obvious invariants 

Ci = ^TtH and C2 = | Tr if^ _ Q Tr (2) 

Note that any normalizable and positive function of Ci and C2 can be regarded as a 
proper distribution function in RMT. To have a specific distribution, one needs to impose 
constraints other than the symmetry consideration. For example, one may impose the 
statistical independence of the variables in H. In the case of Hermitian random matrices, 
this leads to the following probability density function 

P{H) = Afexp {~A Tr H^) , (3) 
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where A/" is a normalization constant. For convenience, we re-parameterize H by using 
Pauli matrices cr and the identity matrix ctq as 

H = eao + X-(T, (4) 

where X = {x,y,z). Parameter e is a trivial parameter because it just reflects the 
reference point of the eigenvalues. In this parametrization, the two U{2) invariants are 

Ci = e and C2 = x'^ + y'^ + . (5) 

A natural U (2)-invariant measure of H is given by 

dfi{H) = dedxdydz. (6) 

A GUE ensemble can then be defined by use of the P{H) in Eq. (I3j) and the d^{H) 
in Eq. ([6]). For the sake of comparison with other extended cases, we fix A = 7r/2 
throughout. Then, after normalizing the total probability to unity, one finds A/" = 7r/2 
as well as the following level-spacing distribution 

where 

s = 2^/C~2 = 2y^x^ + y^ + z^ (8) 

represents the spacing of two eigenvalues. In the three-dimensional parameter space 
spanned by {x,y,z), s can be interpreted as twice of the distance from the origin. 
Because we have fixed A, the average level spacing is not unity here. 

To reduce the symmetry of P{H)dfi{H), one may reduce the symmetry of either 
P{H) or dfi{H). Although in the following we focus on modifications made to dfi{H), 
they can be equally interpreted as modifications to P{H). Consider then 0(2) as 
a reduced symmetry invariant group. Under any 0(2) transformation, there is an 
additional invariant, i.e., 

O3 = y. (9) 

Evidently, in general an arbitrary function / (01,0*2,03) will be invariant under 0(2), 
but not U{2). Hence the following new measure 

dfi'iH) = fiC^,C2,C,)dfiiH) (10) 

will have a reduced symmetry as compared with dfi{H). 

3. Case studies of random Hermitian matrices 

3.1. From GUE to GOE 



To illustrate our methodology, we first show how GOE statistics can naturally emerge 
if we choose /(Oi, 02,03) = ^(Os) = This constraint reduces a Hermitian matrix 
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to a real symmetric form. One then obtains 

J P{H)di2'{H) =N j e-"(^'+"'+^'+"')5(i/)dedxd?/dz 

= U ! e-"("'+^'+"')deda;d2 



AT- / se-^'"l^ds. (Ill 



2.0 

Note that in obtaining the above result we first carry out an integration over the 
delta-constraint and then convert one of the remaining integration variables to s (same 
technique will be used for all the following cases). Using the normalization of total 
probability, we find = 1. Then Eq. ( ITT|) directly yields the standard GOE distribution 
of level spacing, i.e., 

PGOE(s) = |se— (12) 

Interestingly, such a reduction from GUE to GOE is not as trivial as it looks. Because 
the off-diagonal matrix element x — iy can be also parameterized by re"'*^, one might 
naively think that setting the constraint 5(0) will equally reduce GUE to GOE. However, 
such a 5(0) constraint is inappropriate because is not invariant under the pre-specified 
symmetry 0(2). So even with this new parametrization, one should use (5(r sin0) = 5{y) 
to correctly reduce the symmetry from U{2) to 0(2). 

3.2. Planar symmetry 

To demonstrate how our intuitive idea may lead to previously unknown results, we 
consider here the 0(2) invariant function /(Oi, C*2, Oa) = 5(0*3 — Ho) = 5(?/ — yo). In 
the parameter space, this constraint describes a plane perpendicular to the |/-axis. The 
ensemble statistics is then characterized by 

j P{H)di2\H) =U j e--(^'+"'+?^'+^')5(y - yQ)dedxdydz 
= N I e-"(^'+"'+^o+"'Meda;d^ 



poo 

= ^fl / se-^'"/^ds. (13) 

2 J2\yo\ 

Upon fixing the normalization constant Af, we can directly read out a new level-spacing 
distribution function from Eq. ( IT3|) . i.e., 

0, s<2\yo\] 

P,GOEis)={ ^ _ ^ (14) 

s > 2\yo\. 



TT 



-sexp 



-fgGOE(s) can be termed as a "gapped-GOE" distribution because it resembles the GOE 
distribution except that the level spacing must be larger than 2\yo\. The gapped feature 
can be regarded as an extreme situation of level repulsion. Obviously, this result reduces 
to GOE when yo = 0. 
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Certainly, any distribution of y other than 6{y — yo) is still 0(2) -invariant. So we 
may superpose many ensembles with different yo together to form a larger ensemble 
with the same reduced symmetry. Let us consider the averaging of yo over the Gaussian 
distribution P(|/o)d?/o = ee~'^''^^od?/o with e > 0. Then the overall level-spacing statistics 
P{s) is given by a weighted distribution, 

Pis) ^ {P,GOEis))^^ 

/oo 
ee-^^'yoP^ooE{s)dyo 
-oo 

= -ese~^''/^ f'^ e"(i-^')^odyo. (15) 

2 J-s/2 

P{s) obtained above is essentially the well-known interpolation between GUE and GOE 
[5]. This becomes obvious if we set e = 1, resulting 

P{s) = e--^V4 = ^s'e--^y\ (16) 

^ J-s/2 ^ 

which is exactly the Pgue(s) obtained in Eq. ([7]). The case of e = 1 is special because the 
ensemble under this precise condition accidentally possesses a U (2) symmetry invariant 
group. That even the known GOE-GUE distribution can be modeled here further 
indicates that our approach based on the 0(2)-invariant function /(Ci, C2, C3) is simple 
but quite general. 



3.3. Cylindrical Symmetry 

We now start to explore various ensembles of 0(2) invariance by allowing the function 
/(Oi, 02,03) to live on different geometrical objects in the {x,y,z) parameter space. 
Here we let /(Oi, O2, O3) = 6{C2 — 0| — Pg) = S{x'^ + 2"^ ~Po)' which describes a cylinder 
invariant to arbitrary rotation about the y axis. As a result, we have 

J P{H)dij'{H) =U j e-(-'+-'+2^'+-')5(a;2 + - pl)dedxdydz 

= U j e-^^''+p'+y'^{p^ - pDpdpdOdedy (17) 
= A/-7r / e-"^ /Mg, (18) 

where we introduced the polar coordinates: x = p cos 6 and z = p sin 6. 

The normalization condition requires Af = ^e'^'^o. Accordingly the level-spacing 
distribution is found to be 

0, s < 2po; 

» > 2p„. (19) 



exp 



This P{s) in Eq. ([T^ has two noteworthy features. First, it is not only gapped but 
also singular: P{s) diverges as s approaches 2po. Second, as po — )■ 0, P{s) — )■ e~'^^^/^, 
which is a half-Gaussian distribution and hence predicts level clustering. This limit is 
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somewhat expected because as po 0, both x and z will be fixed at zero and hence s 
is solely determined by \y\, whose distribution is half-Gaussian. 

3.4- Parabolic Symmetry 
Here we consider 

f{Ci,C2, Cg) = S{x^ + z^- 2ay) + S{x^ + + 2ay) (20) 

with a > 0. This hence introduces another invariant measure d^'{H). In the parameter 
space, the constraint /(Ci, (72,6*3) describes two paraboloids with the y-axis as their 
symmetry axis. After some calculations we obtain 

J P{H)dfi'iH) = Y j e""f"'+^'^'+"^'^°'lrfMed?/ 

= Afn / -^====e-^y^ds, (21) 
Jo Vs2 + 4a2 

where we still adopt x = p cos 9 and z = p sin 6 as an alternative parameterization. The 
first line of Eq. (121 p indicates that after we confine the ensemble to a parabolic surface, 
the ensemble can still be interpreted as that generated by three statistically independent 
variables e, y, and 9, whose distribution is given by a Gaussian distribution with zero 
mean, an even distribution patched by two Gaussian tails, and a uniform distribution in 
[0, 2n), respectively. In this sense, statistical independence of random variables is again 
respected. 

e~™^ s ^2/4 

^^'^ = Erfc (a0F) V^^tW'^^ ^^^^ 
This distribution represents a smooth interpolation between the GOE distribution 
-Pgoe(s) and a half Gaussian distribution. For a — )■ 00, Pgoe(s) is recovered, whereas 
for a — > 0, we have P{s) — e~'^*^/^. However, it should be also noted that for any fixed 
a 7^ 0, we always have P{s = 0) = and hence level repulsion in the regime of small- 
s. Note also that such an interpolation between level repulsion and level-clustering is 
much different from that of the well-known Berry- Robnik distribution [SlIH]- Based on 
physical considerations of quantum systems with a mixed classical phase space structure, 
the Berry-Robnik distribution interpolates between GOE and the Poisson distribution; 
whereas our result here is within the framework of random matrix ensembles with certain 
0(2) symmetry. In fact, the Poisson distribution never emerges in our calculations. 

Another constraint with a parabolic symmetry can be chosen as /(Ci, (72,6*3) = 
S{x^ + z"^ — ^y'^) with 7 > 0. This constraint describes a surface obtained by rotating 
a parabola about an axis perpendicular to its symmetry axis. The final level-spacing 
distribution function is found to be 

/ P{H)dp'{H) =M r J- e-^''/^ds, (23) 

J Jo ^/^r^^^/^^^^ 

which is somewhat complicated. The P{s) distribution here differs from all previous ones 
in that it always approaches a nonzero value for s — )■ 0, thus predicting level-clustering 
even in the small-s limit. 
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3.5. Conic Symmetry 

Here we let /(Ci, C2, C3) = 5 [x^ + — /?(?/ + yoY] with /3 > 0, which describes a double 
cone in the parameter space. Performing similar calculations as before, we arrive at 

. r ---'-f'-t^°"d.. (24, 

J2g\yo\ - 4^ 

where = + /S). This result is very similar to the cylindrical case discussed 

earlier, so the level-spacing distribution is also both singular and gapped. When yo 
vanishes, we obtain again a half-Gaussian distribution. 



3.6. Discussion 

Our case studies indicate that dramatically different level-repulsion behavior can 
be obtained for Hermitian random matrix ensembles with an orthogonal symmetry 
invariant group. The new distributions found here might be useful in some intermediate 
situations, where some hidden continuous symmetry, albeit smaller than the canonical 
group, plays a crucial role. 

Physical realizations of the ensembles proposed in this section are also possible. 
Let us first imagine an ensemble of spin-| particles embedded in a solid and subject to 
a magnetic field of strength B along the ^/-direction. If we further assume that i? is a 
constant and an internal magnetic field is randomly oriented in the x — z plane, then 
the 2x2 spin Hamiltonians describing this physical situation form an ensemble with a 
planar symmetry. If B is also random, but the internal field strength is fixed, or if there 
exists an induction mechanism that establishes a linear or quadratic dependence of the 
internal field strength upon B, then other types of invariance symmetries discussed in 
this section may be realized. 



4. PT-symmetric random matrix ensembles 

PT-symmetric matrices are non-Hermitian matrices, but they can generate real 
eigenvalues due to the underlying invariance under parity and time reversal operations. 
Recently such systems have attracted a great deal of theoretical and experimental work 
[10] . Here we aim to study how "PT-symmetric random matrix ensembles might be 
constructed and explore the associated level statistics. Because a 2 x 2 PT-symmetric 
matrix has two more parameters than a 2 x 2 Hermitian matrix (detailed below), an 
interesting conjecture is that a 2 x 2 Hermitian random matrix ensemble could be just 
one symmetry-reduced result from a greater symmetry class. If this is true, then more 
types of random matrix ensembles may be identified. 

We start with a generic 2x2 non-Hermitian matrix characterized by eight real 
parameters, 

rr _ ( an + ibn au + ibu 

-n general I . . 

\ 021 + IO2I ^22 + IO22 
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Consider now similarity transformations, i.e., GL{2,C) as the canonical group. Under 
a similarity transformation, the invariants are 

Ci = |TriJ and C2 = | Tr //^ - (| Tri^)^ (26) 

Note that both Ci and C2 are complex quantities. The real parts and the imaginary 
parts are independently invariant under GL{2, C). The intuitive measure d/x(ifgenerai) = 
n?j=i dajjdfejj is also invariant under GL{2,C). 

To get a specific distribution function for the matrix elements, we need to impose 
constraints on top of the symmetry consideration. In line with the GUE and GOE cases, 
here we start from a statistics-independence consideration (which however depends 
on the parametrization of -^general)- Under our parameterization detailed above, 
if we request that -P(-f^generai)d/i(-ffgenerai) describcs the parameters in a statistically 
independent manner, then P(-ffgenerai) is found to be trivial because it depends only on Gi 
(after some lengthy calculations). For this reason we introduce a new parameterization 
by expanding i^generai in Pauli matrices, i.e., 

i^general = (c + ie)(To + (X + iR) ■ CT, (27) 

where X = (x, y, z) and R = (p, g, r). With this set of new variables, the above GL{2, C) 
invariants become 

Ci = e + ie and = X^ - R^ + 2iX ■ R, (28) 

and the invariant measure becomes d/i(ifgenerai) = dededXdR. The requirement that 
this set of new variables should be statistically independent finally leads to 

P(i^generai) = AfP{e)P{e) cxp [-2A (X^ - R^)] , (29) 

where -P(e) and -P(e) are arbitrary distribution functions and A is a constant. As a 
convenient choice, we choose -P(e) and -P(e) such that 

P ( Hgenersd ) d/i ( -ff general ) 

= ATexp [-ARe (Tr^r|„i)] dededXdR 

= Afexp [-2 A (e^ - e2 + - R^)] dededXdR, (30) 

where is a normalization constant where appropriate. Such type of 

P(ifgenerai)d/i(-ffgenerai) appears uaturally from our consideration of statistical indepen- 
dence. At this point 

P(i7generai)d/i(-ffgenerai) is much different from what is widely used in previous studies 
of non-Hermitian random matrix ensembles [21 [TJ [11], |12]. In terms of the above 
parametrization, the previously considered distribution function [21 [TJ HH [12] has the 
form 

-P(-^^non— Hermitian) A/'cxp A Tr ^-ffjjQjj_j£gj.jjjj^^jg^jj-ffnon— Hermitian^ 

= Arexp[-2yl(e2 + e2 + x2 + R2)] (31) 
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Our P(-ffgenerai) IS invariant upon GL{2, C) transformations and the -P(-ffnon-Hermitian) 
in Eq. fl3Tl) is only U{2) invariant. To verify this, it is easy to show that the individual 
length of X or R is invariant if only U (2) transformations are considered: 

C3 = X2 and C4 = R^. (32) 

Therefore, the -P(-ffnon-Hermitian) in Eq. (I3T!) can be considered as a result of symmetry 
reduction by applying the U (2) symmetric constraint 

/(Ci, ^2, C3, C4) = exp [-AA (e^ + R^)] 

to Eq. fl30l) . Certainly, as we will see below, this is only one special example U{2) 
symmetric constraint. 

The probability density -P(-ffgenerai) and the invariant measure d/i(ifgenerai) we 
proposed above still do not suffice to define an acceptable random matrix ensemble 
with real eigenvalues. For all eigenvalues to be real, a 2 x 2 matrix must have the VT 
symmetry |6]. This is the case under the following conditions: 

ImCi = e = and ImC2 = 2X-R = 0. (33) 

Note that these constraints are GL(2, C)-invariant also, so we are not reducing the 
invariance symmetry yet. Following exactly the same line as in the Hermitian cases, we 
implement these two constraints by modifying the invariant measure to the following: 

d/i(ifpr) = S{e)6 j dededXdR. (34) 

The denominator in the second (5-function is also GL{2, C)-invariant and here it is to 
give a proper dimension. 

Equations (!30|) and now define an ensemble of 2 x 2 PT-symmetric random 
matrices. Purely for the sake of further calculations, we now adopt the parametrization 
similar to the one used in Ref. |6] for PT-symmetric matrices, namely, 

H-pj- = eo-Q + (th*" + iu sin rjn^ + iu cos rj n"^) ■ cr, (35) 

with 

n^' = (sin 6 cos if, sin 6 sin ip, cos 6) , 

= (cos 9 cos (f, cos 9 sin ip, — sin 9) , 
n*^ = (— siny9, cos 0). (36) 

In term of the new variables, x = 7 sin 6' cosy?; y = 7 sin ^ sin y?; z = ycos9; p = 
u {sin rj cos 9 cos ip — cos rj sin ip); q = i/ (sin 77 cos 6* sin + cosrj cos ip); t = — z/sinr/sin^^, 
with X^ = 7^ and R^ = i^^. With this parameterization, we have 



P{HvT)d^iiH-pT) 

M j exp [-2 A (e^ + 7^ - v"^)] 7z/sin^^7^^^ded7dz/dMy?dr7. (37) 
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To focus on real eigenvalues, we next introduce a constraint such that only VT- 
symmetric matrices with real eigenvalues are included, which can be realized by the 
constraint 7 > [6]. We hence introduce a GL(2, C)-invariant constraint as a step 
function 9, i.e.. 



With this rather formal expression, a simple analysis indicates that this is still not a 
normalizable probability distribution (A/" — t- 0). The apparent reason is that the factor 
exp [—2A (e^ + 7^ — u"^)] in PlH-pf) does not decay quickly with increasing 7^ so long 
as + 7^ — z/^ is kept small. A deeper explanation should be linked with the non- 
compactness of the GL{2, C) group, which is closely related to the Lorentz group 0(3, 1). 
That is, to have a normalizable distribution function, we need to reduce the invariance 
symmetry to a compact subgroup. In this case, the largest compact subgroup is U{2). 
This understanding is also confirmed below when we introduce constraints to reduce the 
invariance symmetry to U{2) (note: this does not mean the matrices themselves have 
to be complex Hermitian!). 

4.I. From VT symmetry to Hermiticity 

Obtaining Hermitian matrices from PT-symmetric matrices can be done by a constraint 
= z/^ = 0. This constraint is U{2) invariant because it is about the square of the 
length of R. Imposing this constraint, we obtain the following f/(2)-invariant ensemble 



where the ?7-parameter becomes irrelevant and can hence be integrated out. Since here 
the level spacing s is just 27, the above expression immediately leads to the standard 
GUE distribution Pgue(s) given in Eq. ([7]). Indeed, using our early variables for 
Hermitian matrices, the expression for J P{HTy'j-)dfi{H-pf) is also seen to be identical 
with f P(-ffHermitian)d/i(-ffHemiitian)- Thus, a Standard Gaussiau unitary ensemble of 2 x 2 
Hermitian random matrices is seen to be a special example of our "PT-symmetric random 
matrix ensemble under a U{2) constraint. 

4-2. General constraint on v 



M / exp {-2 A {(? + 7^ - z/2)] 7z/sin^v^72 - v'^Q - 




) 



X ded7di^d^d<^d?7. 



(38) 




(39) 



Our success in reducing P {H'pr)(^n{H-pr) to P(-ffHermitian)d/u(ifHermitian) suggests the 
existence of [/(2)-invariant ensembles not seen before by imposing more general U{2) 
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constraints. As a very simple example, we may consider a slice of v — uq ^ of our 
original "PT-symmetric ensemble to generate an interesting sub-ensemble with again 
U{2) symmetry. With this new constraint we have 

J P{Hrr)MHvT) 

^Af J exp [-2A (e^ + 7^ - u^)] jiysm0y/j^ - u'^Q{j^ - iy^)S{iy - uo) 
X ded7di/d6'd</7dr?. (40) 

After carrying out the necessary integrals, one finds that the statistics of the level 
spacing s = 2^/7^ — Uq is again given by Pgue(s), i-e., the standard GUE result. This is 
somewhat unexpected because after all, our ensemble consists of non-Hermitian random 
matrices! 



4.3. Constraint on 7 

In exactly the same fashion, we may introduce a constraint about 7 to generate new 
ensembles. Because 7 itself is again C/(2)-invariant, an arbitrary function of 7 may be 
used as a constraint, producing many types of [/(2)-invariant ensembles in a manner 
analogous to our previous Hermitian ensembles. Here we consider the simplest case, in 
which the constraint is given by 7 = 70. Then, we can construct the following ensemble, 

J P{Hvt)MHvt) 

= exp [-2 A (e^ + 7^ - v'')] 7i/sin^V7^^^©(7' - - 7o) 

X ded7di/d^d</7d77. (41) 
Integrating out all the free variables under a fixed level spacing s — 2\/7o — v"^, we find 

WW = I '^ ^'"^ (42) 

1 0' ^ > 270- 

This level-spacing statistics can be called as a "truncated-GUE" distribution because it 
is exactly the same as the GUE distribution for s less than a threshold value and zero 
elsewhere. With this result, it becomes apparent that many interesting level-spacing 
statistics may emerge if we consider further a distribution of 7 values. This will not be 
pursued here. 



5. Conclusion 



In conclusion, by reducing the invariance symmetry of random matrix ensembles via 
various constraint functions, our straightforward calculations are able to reveal novel 
features in the level-spacing statistics of 2 x 2 random matrices, including both Hermitian 
matrices and PT-symmetric non-Hermitian matrices. In the Hermitian case, we find 
singular distributions, half Gaussian distribution, "gapped-GOE" distribution, as well 
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as distributions interpolating between GOE and half-Gaussian distributions. Our 
construction of random 2x2 Hermitian matrix ensembles is remarkably different from 
a recent interesting study by Berry and Shukla [13], where intriguing level-spacing 
distribution functions for generalized Gaussian ensembles of 2 x 2 real symmetric random 
matrices are found in the absence of 0(2) invariance symmetry. In the non- Hermitian 
case, our symmetry-reduction approach to PT-symmetric random matrix ensembles 
appears to be more general than previous treatments of non-Hermitian random matrix 
ensembles [21 [TJ HH [12] (which can be also regarded as an outcome of symmetry 
reduction), thereby establishing new types of 2 x 2 non-Hermitian random matrix 
ensembles with ?7(2)-invariance and unexpected level-spacing statistics. 

It is an open question whether results based on 2 x 2 matrices are relevant to VT- 
symmetric random matrices with a very high dimension. We do not speculate on this 
challenging issue because at present how to parameterize such matrices is still unclear. 
What can be anticipated, however, is that many other types of random matrix ensembles 
can be considered once the parameterization of high-dimensional PT-symmetric random 
matrices are known, with the known GOE and GUE cases as a very small subset of a 
much broader class of random matrix ensembles. 
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